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Using an exact diagonalization technique within a generahzed Mott-Hubbard Hamiltonian, we 
predict the existence of a ground state persistent current in coherent two-dimensional semiconductor 
quantum dot arrays pierced by an external magnetic flux. The calculated persistent current, which 
arises from the nontrivial dependence of the ground state energy on the external flux, exists in 
isolated arrays without any periodic boundary condition. The sensitivity of the calculated persistent 
current to interaction and disorder is shown to reflect the intricacies of various Anderson-Mott- 
Hubbard quantum phase transitions in two dimensional systems. 



X3 
O 

o 



> 
o 

(N 
O 

On 
-I— > 



I 

O 

o 



X 
S3 



Significant recent advances in semiconductor materi- 
als growth and nanolithography techniques now allow 
101 for the fabrication of multiple quantum dot systems 
in modulation-doped two dimensional electron struc- 
tures which are quantum mechanically coherent, and 
may therefore be considered two dimensional "artifi- 
cial molecules" (with individual quantum dots being the 
"atomic" constituents of this artificial molecule). The- 
oretical work on multi-dot systems has mostly concen- 
trated on the two limiting situations (coherent dots 
with no Coulomb interaction and Coulomb blockade 
of individual dots |^). In this Letter we consider coher- 
ent two dimensional (2D) arrays of quantum dots (i.e. 
2D artificial molecules) taking into account jJJ^ both 
quantum fluctuations arising from interdot hopping and 
electron-electron interaction effects through a generalized 
Mott-Hubbard Hamiltonian Our main result is the 
prediction of an equilibrium ■persistent current in finite 
2D arrays of semiconductor quantum dots in the pres- 
ence of an applied magnetic field transverse to the 2D 
plane. A measurement of this persistent current in small 
2D arrays will be definitive evidence for the formation of 
an artificial quantum dot molecule. The predicted per- 
sistent current is a periodic function of the external flux 
and exhibits clear signature of the competition between 
coherence and interaction in artificial 2D molecules. We 
include disorder effects in our calculation and show that 
our predicted equilibrium persistent current should be ex- 
perimentally observable in currently (or soon-to-be) ac- 
cessible 2D quantum dot arrays. 

It is well-known j^-^ that gauge invariance and single- 
valuedness of electron wavefunction allow for the exis- 
tence of a ground state persistent current in normal metal 
rings surrounding an external magnetic flux. The intrin- 
sic magnetic moment associated with this persistent cur- 
rent, which is proportional to the current in one dimen- 
sional (ID) rings, is an oscillatory function of the external 
flux with a period equal to the elementary flux quantum 

= h/e. The existence of such an oscillatory persistent 
current in normal metal rings has been experimentally 



verified |9|-|l 1[] . Our predicted ground state persistent 
current in 2D quantum dot arrays has some significant 
differences with the existing persistent current theoret- 
ical analyses in the ID ring geometry. The system we 
consider is a finite coherent 2D quantum dot array where 
the one electron states are in the atomic tight binding 
limit (i.e. the electron states tend to be strongly local- 
ized on individual quantum dots except for weak inter- 
dot quantum tunneling) in contrast to metallic systems 
which are nearly free electron like. In addition, we con- 
sider a finite 2D array without any periodic boundary 
condition in contrast to the ID ring geometry, which, by 
definition, is periodic. Another significant aspect of our 
prediction is the inclusion of quantum tunneling, disor- 
der, and Coulomb interaction effects on equal footings 
because in real semiconductor quantum dot arrays all of 
these effects are expected to be important. 

We model the finite square 2D quantum dot array of 
size Lx X Ly = L dots through the following physically 
motivated generalized 2D Mott-Hubbard type strongly 
correlated Hamiltonian Q (c, are electron field oper- 
ators and n is the density operator): 
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where 



bij = I A ■ lij is the usual Peierls gauge phase 
factor arising from electron hopping on a lattice in a 
transverse external magnetic flux (with A the vector po- 
tential and lij the spatial vector along a plaquette) . The 
indices i, j are the quantum dot spatial positions in the 
array and a denotes an electron state on the quantum 
dot. The generalized Mott-Hubbard parameters U, V 
and t have the usual significance of the on-site Coulomb 
interaction, the long range part of the Coulomb inter- 
action and the hopping kinetic energy term respectively 
within the tight binding description of the 2D array. We 
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consider two spin-split energy levels per dot due to 
quantum confinement in a dot. We assume the magnetic 
field to be weak enough for single-particle level crossing 
effects not to be important. The Coulomb interaction in 
the array is expressed in terms of the capacitance ma- 
trix representation ^ . The quantum dot array is in the 
tight binding atomic limit, and we therefore keep only 
the nearest neighbor tunneling term in tij = t. Within 
the model defined by Eq. (1), U, V, t are parameters 
to be given as inputs using a physical description of the 
quantum dot array. 

We are interested in calculating the equilibrium persis- 
tent current flowing in the array at T = in the presence 
of an external flux threading the 2D lattice. It is easy to 
show by using commutation properties of H with the po- 
larization operator that the transverse component of the 
magnetic moment density operator, m^, is given by the 
gauge invariant result 



where $ is the total uniform magnetic flux piercing the 
finite 2D system. (We note that the external magnetic 
field is in the z direction and with the 2D array in the x—y 
plane, and therefore is the only non-zero component 
of the magnetization density.) A non-zero value of the 
static magnetization density, as measured by a finite , 
is equivalent to the existence of a ground state persistent 
current (whose strength is exactly the same as the mea- 
sured value of rriz because of the two dimensional nature 
of the problem) . Our goal is to calculate the equilibrium 
magnetic moment density for 2D quantum dot arrays us- 
ing Eqs. (1) and (2). An experimental measurement of a 
finite m,z (as a function of external fiux) will be a direct 
verification of the quantum dot molecule formation. 

To obtain the equilibrium persistent current, / = m^^ 
we calculate the ground state energy by doing an exact 
diagonalization of H, defined by Eq. (1), in the basis of 
the total number of electrons N in the array and the to- 
tal spin Sz in the direction of the field, and then perform 
minimization over Sz to find the stable ground state. The 
magnetic moment density in a finite 2D quantum dot ar- 
ray is calculated using Eq. (2). Existing semiconductor 
nanofabrication techniques are at best capable of produc- 
ing coherent arrays which are 2 x 2 to 3 x 3 quantum 
dots in size. The Lanczos diagonalization technique that 
we use is capable of giving the exact interacting ground 
state for upto a 3 x 3 array of our interest. We calcu- 
late the persistent current in interacting 2D arrays upto 
i = 3 X 3 sizes for all band fillings n{= N/2L) by vary- 
ing the number of electrons N in each system between 
N = 1 and N = L. We use realistic Mott-Hubbard pa- 
rameters which approximately correspond to 2D GaAs 
quantum dot systems (at T = 0). We choose the 
intradot single particle level spacing A = 0.3 meV to 
be comparable in magnitude to the interdot hopping pa- 
rameter t = 0.1 meV, with the intradot Coulomb charg- 
ing energy U ^ I meV. With this choice of parame- 



ters the 2D quantum dot array is in the interesting co- 
herent "molecular" state referred to as the "Collective 
Coulomb Blockade" regime , where both coherence and 
interaction play comparable non-perturbative roles, and 
the Coulomb blockade of individual dots is destroyed. 
Disorder is included in our calculation through a spin- 
independent parameter W ^ A which denotes the half- 
width of a uniform distribution of random on-site quan- 
tum dot energies centered around A. The spin splitting 
(set to be 0.03A) of the single particle energy levels is 
assumed to be independent of W. 
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FIG. 1. Typical calculated persistent current for three dif- 
ferent values of band filling (n = l/6(a), 5/18(6), l/2(c)) in 
the 3x3 array is plotted versus the magnetic flux through the 
elementary cell in units of (j>o for six different sets of U, V, W 
(see the legend in the Fig. 2 for the key to each line in the 
plot) in the 3x3 system. 

We show our typical calculated persistent current as a 
function of the external flux in Fig. 1 for three different 
values of band filling (n = 1/6, 5/18, 1/2) in the 3 x 3 
system. Results for six different situations (as described 
in the figure captions) for various values U, V, and W are 
shown for each filling with t, A fixed throughout the cal- 
culation. The flux (j) is through each elementary cell (or 
plaquette) of the 3x3 system, and the calculated equi- 
librium persistent current is in absolute units. The most 
signiflcant feature of Fig. 1 is the existence of a persistent 
current, which is periodic in 4>/4'o {note that the total flux 
through the system is $ = 40 by virtue of the four equiv- 
alent plaquettes in the 3x3 quantum dot array), whose 
magnitude is reasonably large (> nA). The magnitude 
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of the current depends in a rather comphcated way on 
the Mott-Hubbard parameters V and W , and also on 
the band-filhng (in particular, on whether n = 1/2 or 
not). We can, however, make the following qualitative 
observations: (1) Away from half- filling, finite disorder 
is the most prominent destructive effect on the current 
amplitude, with finite interaction ([/, V ^ Q) producing 
some enhancement of the disordered current but never to 
the free electron (J7, W — Q) value. Without any 
disorder {W = 0), finite U and V have reasonably small 
effects on the current magnitude away from half-filling. 
(2) At half-filling, however, finite on-site Coulomb repul- 
sion ([/ ^ 0) by itself dramatically suppresses the per- 
sistent current with the long-range Coulomb interaction 
{V 7^ 0) opposing this dramatic suppression due to finite 
U. (3) Disorder, by itself {W ^ Q, U, V =0), seems 
to produce similar quantitative suppression of the per- 
sistent current at all fillings except at half-filling where 
disorder enhances (slightly) the interacting clean system 
current. (4) Long range Coulomb repulsion {V ^ 0) al- 
ways enhances the current magnitude from its finite U 
and W values. (5) The various structures (jumps and 
discontinuities) in the current arise from many-body en- 
ergy level crossings which occur in the system at specific 
fillings and values of the magnetic flux. The discontinu- 
ities of the persistent current seen in Fig. f correspond 
to the change of the total orbital momentum in the 2D 
array, which can be accompanied by spin-fiip transitions 
in the system. Finite disorder {W ^ 0) smoothens the 
discontinuities in the noninteracting [U, V = Q) current, 
and removes some of the discontinuities in the interact- 
ing (J7, V ^ Q) current. There are no discontinuities 
or spin flips at half-flUing in the interacting system (Fig. 
1(c)). (6) The sign of the persistent current at zero mag- 
netic flux is generally dependent on the electron filling 
and the geometry of the array. We find the interacting 
{U ^Q, V, W = 0) persistent current in 2 x 2, 3 x 2, 
4x2, and 3x3 arrays to be diamagnetic for n — 1/2L 
and paramagnetic for n — 1/2. The random disorder of 
strength W (= A, as used in our calculations) does not 
change the sign of the persistent current. (7) The current 
distribution in the clean 2D array is determined by the 
spatial symmetries of the system and band filling in the 
array. In the clean 3x3 lattice the current flows along the 
boundary for all electron fillings in the system. However, 
the persistent current distribution is disorder realization 
dependent. 

The most important feature of our results, the strong 
suppression of the persistent current at half-filling by the 
on-site Coulomb repulsion [/, is a direct finite size mani- 
festation of the Mott-Hubbard metal-insulator transition, 
which is a true phase transition in the thermodynamic 
limit. The long-range Coulomb repulsion {V ^ 0) op- 
poses this transition in a finite system (leading to an 
enhancement in the finite {/-suppressed current ampli- 
tude for V y^Oat half-filling). Effects of finite V and W 
on the Mott-Hubbard transition in a 2D Hubbard model 
(in the thermodynamic limit) arc not rigorously known. 



but our results are consistent with the expectation that, 
in general, finite values of V and W should oppose the 
metal-insulator Mott transition at half-filling. To further 
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FIG. 2. A plot of the RMS persistent current {J^}^/^ 
versus the number of electrons A*' in the 3x3 array for the 
six different sets of U, V, W (as shown in the legends) in the 
3x3 system. AH energies are given in the meV units. 

quantify the relative magnitudes of the effects of finite 
U, V, W on the persistent current at various band fill- 
ings, we show in Fig. 2 our calculated root mean square 
current amplitude as a function of the number of elec- 
trons in the 3x3 quantum dot array for various values of 
Mott-Hubbard parameters. The rms current magnitude, 
(/2^i/2^ is obtained by averaging the typical current (Fig. 
I) over one flux period. The six different parameter sets 
with different values of U, V, W in Fig. 2 clearly show 
the dominant effect of the Mott-Hubbard transition (for 
finite U) at half-filling (TV = 9), the generic destructive 
effect of finite disorder at all fillings, and the generic ten- 
dency of the long-range Coulomb interaction homogenize 
the electron density in the system and consequently 
to enhance the persistent current. 

In the well-studied single ID ring geometry, the effect 
of disorder (in the absence of any Coulomb interaction) 
on the persistent current is an exponentially strong sup- 
pression [T^ of the current arising from the Anderson 
localization phenomenon associated with all electronic 
states being exponentially localized in one dimension in 
the presence of any finite disorder. For our 2D quantum 
dot array the effect of finite disorder {W ^ 0) is subtle 
(and much softer than in one dimension) because local- 
ization effects are logarithmically weak in two dimen- 
sions for small W. (In fact, in the presence of an external 
magnetic field breaking the time reversal symmetry, the 
localization effect is even weaker |l4| and may be hard 
to discern in a finite sample.) This explains the rather 
benign effect of having a finite W in Figs. 1 and 2 of our 
results, which should be contrasted with ID ring calcula- 
tions where the persistent current is strongly suppressed 
by finite values of W. To better quantify our finite disor- 
der results we show in Fig. 3 a log-log plot at half-filling 
of the rms current {Py^^ without any interaction effects 
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FIG. 3. A log-log plot of {Py'^ versus W averaged 
over 100 disorder realizations for each value of W is shown 
(f/ = 1/ = 0) at half-filling in 3 x 3 (crosses), 4x4 (tri- 
angles), 5x5 (solid line), 6x6 (pluses). Asterisks show the 
interacting {U = 1 meV, V = 0) results at half-filling for the 
3x3 system. Inset: (J^)^''^ at half-filling versus system sizes 
(L): triangles (2D array), squares (ID ring). (The average is 
taken over the range of the same total flux through the sys- 
tems, which varies from —(f)o/2 to 00/2.) The dashed lines in 
the plot are the best linear fits to the corresponding data as 
described in the text. 

(t/, y = 0), averaged over 100 disorder realizations for 
each value of W, as a function of the disorder strength 
W for various system sizes (3x3, 4x4, 5x5, 6x6). In 
plotting these results, we have factored out a scale factor, 
ric = (-L^/^ — \Y /L, so that the results for various sys- 
tem sizes fall on top of each other, showing approximate 
current scaling with system size and disorder. (The scale 
factor Tic = {L^/"^ — l)^/i is the number of unit cells 
or plaquettes in each square array of size L.) All four 
non-interacting disordered results scale very well whereas 
the one disordered result with interaction deviates from 
the universal scaling at low disorder strength where the 
Mott metal-insulator transition effect dominates. The 
two dashed straight lines in Fig. 3 give the best fits to 
weak and strong disorder scaled currents, leading to the 
following empirical results for the effect of disorder on 
the persistent current: {P)^ ~ {L/nc) g(W^), with the 
scaling function g{W) being given by 

(W\ / W^"'', 7 = (6-4±2.8) X 10-2, T4^< 1.55 7rt 

The very small value of the scaling exponent 7 in the 
weak disorder limit is consistent with the expected loga- 
rithmic weak localization ijl^ in a 2D array, which even- 
tually crosses over to strong localization for large W . 
We also expect interaction effects to be particularly im- 
portant at weak disorder, but not so at high disorder, 
which is what is seen in Fig. 3 even at half-filling. Weak 
disorder, in fact, produces anti-localization effect (Fig. 
1(c) and 2) at half-filling in the presence of finite U, by 
slightly enhancing the persistent current from its finite U - 



suppressed value. Whether these last two findings remain 
true in the thermodynamic limit is, however, unknown. 

Note that the scale factor Uc in Fig. 3 indicates that 
the persistent current scales with the perimeter size of 
the array or the square root of the number of dots, L^^^, 
in the system. We have explicitly verified this system size 
dependence by doing calculations for various system sizes 
at various fillings. As an example, we show as an inset 
of Fig. 3 a comparison of the system size dependence 
of the rms current in ID ring and 2D array geometries, 
finding an approximately linear dependence on L" with 
the best fit a « 0.46 (1.1) in the 2D array (ID ring). 
Thus, the persistent current is proportional to the length 
of the boundary in both ID ring and 2D array geometry. 

In conclusion, we predict, based on an exact diagonal- 
izaton study of the generalized disordered Mott-Hubbard 
Hamiltonian, the existence of a ground state persistent 
current in finite 2D quantum dot arrays which should be 
experimentally observable. Direct observation of an equi- 
librium persistent current in coherent 2D quantum dot 
arrays will not only verify the formation of an artificial 
quantum dot molecule but will also shed light on the in- 
terplay among coherence, disorder, long- and short- range 
Coulomb interaction effects on Anderson-Mott-Hubbard 
quantum phase transitions. 

The authors thank Charles Stafford for useful discus- 
sions during the early part of this work. This work was 
supported by the U.S. - O.N.R. 
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Using an exact diagonalization technique within a generalized Mott-Hubbard Hamiltonian, we 
predict the existence of a ground state persistent current in coherent two-dimensional semiconductor 
quantum dot arrays pierced by an external magnetic flux. The calculated persistent current, which 
arises from the nontrivial dependence of the ground state energy on the external flux, exists in 
isolated arrays without any periodic boundary condition. The sensitivity of the calculated persistent 
current to interaction and disorder is shown to reflect the intricacies of various Anderson-Mott- 
Hubbard quantum phase transitions in two dimensional systems. 



Significant recent advances in semiconductor materi- 
als growth and nanolithography techniques now allow 
[1] for the fabrication of multiple quantum dot systems 
in modulation-doped two dimensional electron struc- 
tures which are quantum mechanically coherent, and 
may therefore be considered two dimensional "artifi- 
cial molecules" (with individual quantum dots being the 
"atomic" constituents of this artificial molecule). The- 
oretical work on multi-dot systems has mostly concen- 
trated on the two limiting situations [2,3] (coherent dots 
with no Coulomb interaction [2] and Coulomb blockade 
of individual dots [3]). In this Letter we consider coher- 
ent two dimensional (2D) arrays of quantum dots (i.e. 
2D artificial molecules) taking into account [4,5] both 
quantum fiuctuations arising from interdot hopping and 
electron-electron interaction effects through a generalized 
Mott-Hubbard Hamiltonian [4]. Our main result is the 
prediction of an eqmUhnum persistent current in finite 
2D arrays of semiconductor quantum dots in the pres- 
ence of an applied magnetic field transverse to the 2D 
plane. A measurement of this persistent current in small 
2D arrays will be definitive evidence for the formation of 
an artificial quantum dot molecule. The predicted per- 
sistent current is a periodic function of the external fiux 
and exhibits clear signature of the competition between 
coherence and interaction in artificial 2D molecules. We 
include disorder effects in our calculation and show that 
our predicted equilibrium persistent current should be ex- 
perimentally observable in currently (or soon-to-be) ac- 
cessible 2D quantum dot arrays. 

It is well-known [6-8] that gauge invariance and single- 
valuedness of electron wavefunction allow for the exis- 
tence of a ground state persistent current in normal metal 
rings surrounding an external magnetic fiux. The intrin- 
sic magnetic moment associated with this persistent cur- 
rent, which is proportional to the current in one dimen- 
sional (ID) rings, is an oscillatory function of the external 
fiux with a period equal to the elementary fiux quantum 
= h/e. The existence of such an oscillatory persistent 
current in normal metal rings has been experimentally 



verified [9-11]. Our predicted ground state persistent 
current in 2D quantum dot arrays has some significant 
differences with the existing persistent current theoret- 
ical analyses in the ID ring geometry. The system we 
consider is a finite coherent 2D quantum dot array where 
the one electron states are in the atomic tight binding 
limit (i.e. the electron states tend to be strongly local- 
ized on individual quantum dots except for weak inter- 
dot quantum tunneling) in contrast to metallic systems 
which are nearly free electron like. In addition, we con- 
sider a finite 2D array without any periodic boundary 
condition in contrast to the ID ring geometry, which, by 
definition, is periodic. Another significant aspect of our 
prediction is the inclusion of quantum tunneling, disor- 
der, and Coulomb interaction effects on equal footings 
because in real semiconductor quantum dot arrays all of 
these effects are expected to be important. 

We model the finite square 2D quantum dot array of 
size Lx X Ly = L dots through the following physically 
motivated generalized 2D Mott-Hubbard type strongly 
correlated Hamiltonian [4] (c, are electron field oper- 
ators and h is the density operator): 

i ,a i ij 

^ ^ a ,ij ^ ia^j a ~1~ /i.e.), (1) 

<i,j>,a 

where (j>ij = |- f.j A ■ is the usual Peierls gauge phase 
factor arising from electron hopping on a lattice in a 
transverse external magnetic fiux (with A the vector po- 
tential and lij the spatial vector along a plaquette). The 
indices i, j are the quantum dot spatial positions in the 
array and a denotes an electron state on the quantum 
dot. The generalized Mott-Hubbard parameters U, V 
and t have the usual significance of the on-site Coulomb 
interaction, the long range part of the Coulomb inter- 
action and the hopping kinetic energy term respectively 
within the tight binding description of the 2D array. We 
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consider two spin-split energy levels £„ per dot due to 
quantum confinement in a dot. We assume the magnetic 
field to be weak enough for single-particle level crossing 
effects not to be important. The Coulomb interaction in 
the array is expressed in terms of the capacitance ma- 
trix representation [3]. The quantum dot array is in the 
tight binding atomic limit, and we therefore keep only 
the nearest neighbor tunneling term in tij = t. Within 
the model defined by Eq. (1), U, V, t are parameters 
to be given as inputs using a physical description of the 
quantum dot array. 

We are interested in calculating the equilibrium persis- 
tent current fiowing in the array at T = in the presence 
of an external fiux threading the 2D lattice. It is easy to 
show by using commutation properties of H with the po- 
larization operator that the transverse component of the 
magnetic moment density operator, , is given by the 
gauge invariant result 



dH 
"5$"' 



(2) 



where $ is the total uniform magnetic fiux piercing the 
finite 2D system. (We note that the external magnetic 
field is in the z direction and with the 2D array in the x—y 
plane, and therefore is the only non-zero component 
of the magnetization density.) A non-zero value of the 
static magnetization density, as measured by a finite , 
is equivalent to the existence of a ground state persistent 
current (whose strength is exactly the same as the mea- 
sured value of because of the two dimensional nature 
of the problem). Our goal is to calculate the equilibrium 
magnetic moment density for 2D quantum dot arrays us- 
ing Eqs. (1) and (2). An experimental measurement of a 
finite (as a function of external fiux) will be a direct 
verification of the quantum dot molecule formation. 

To obtain the equilibrium persistent current, I = ni;,, 
we calculate the ground state energy by doing an exact 
dtagonaUzatton of H , defined by Eq. (1), in the basis of 
the total number of electrons N in the array and the to- 
tal spin Sz in the direction of the field, and then perform 
minimization over to find the stable ground state. The 
magnetic moment density in a finite 2D quantum dot ar- 
ray is calculated using Eq. (2). Existing semiconductor 
nanofabrication techniques are at best capable of produc- 
ing coherent arrays which are 2 x 2 to 3 x 3 quantum 
dots in size. The Lanczos diagonalization technique that 
we use is capable of giving the exact interacting ground 
state for upto a 3 x 3 array of our interest. We calcu- 
late the persistent current in interacting 2D arrays upto 
_L = 3 X 3 sizes for all band fillings n{= N/2L) by vary- 
ing the number of electrons N in each system between 
N = I and N = L. We use realistic Mott-Hubbard pa- 
rameters which approximately correspond to 2D GaAs 
quantum dot systems (at T = 0). We choose [4] the 
intradot single particle level spacing A = 0.3 meV to 
be comparable in magnitude to the interdot hopping pa- 
rameter t = 0.1 meV , with the intradot Coulomb charg- 
ing energy ~ 1 meV . With this choice of parame- 



ters the 2D quantum dot array is in the interesting co- 
herent "molecular" state referred to as the "Collective 
Coulomb Blockade" regime [4], where both coherence and 
interaction play comparable non-perturbative roles, and 
the Coulomb blockade of individual dots is destroyed. 
Disorder is included in our calculation through a spin- 
independent parameter VF ~ A which denotes the half- 
width of a uniform distribution of random on-site quan- 
tum dot energies centered around A. The spin splitting 
(set to be 0.03A) of the single particle energy levels is 
assumed to be independent of W . 
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FIG. 1. Typical calculated persistent current for three dif- 
ferent values of band fillingi (n = l/6(a), 5/18(6), l/2(c)) in 
the 3x3 array is plotted versus the magnetic flux through the 
elementary cell in units of i/iq for six different sets of U, V, W 
(see the legend in the Fig. 2 for the key to each line in the 
plot) in the 3x3 system. 

We show our typical calculated persistent current as a 
function of the external fiux in Fig. 1 for three different 
values of band filling (n = 1/6, 5/18, 1/2) in the 3x3 
system. Results for six different situations (as described 
in the figure captions) for various values U, V, and W are 
shown for each filling with t, A fixed throughout the cal- 
culation. The fiux (f) is through each elementary cell (or 
plaquette) of the 3x3 system, and the calculated equi- 
librium persistent current is in absolute units. The most 
significant feature of Fig. 1 is the existence of a persistent 
current, which is periodic in (fi/ctio (note that the total fiux 
through the system is $ = 4(/) by virtue of the four equiv- 
alent plaquettes in the 3x3 quantum dot array), whose 
magnitude is reasonably large (> nA). The magnitude 
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of the current depends in a rather complicated way on 
the Mott-Hubbard parameters U, V and W , and also on 
the band-filling (in particular, on whether n = 1/2 or 
not). We can, however, make the following qualitative 
observations: (1) Away from half-filling, finite disorder 
is the most prominent destructive effect on the current 
amplitude, with finite interaction (U, V ^ 0) producing 
some enhancement of the disordered current but never to 
the free electron (U, V, W = 0) value. Without any 
disorder (W = 0), finite U and V have reasonably small 
effects on the current magnitude away from half-filling. 
(2) At half-filling, however, finite on-site Coulomb repul- 
sion (U 7^ 0) by itself dramatically suppresses the per- 
sistent current with the long-range Coulomb interaction 
(V 7^ 0) opposing this dramatic suppression due to finite 
U. (3) Disorder, by itself (W ^ 0, U, V = 0), seems 
to produce similar quantitative suppression of the per- 
sistent current at all fillings except at half-filling where 
disorder enhances (slightly) the interacting clean system 
current. (4) Long range Coulomb repulsion (V ^ 0) al- 
ways enhances the current magnitude from its finite U 
and W values. (5) The various structures (jumps and 
discontinuities) in the current arise from many-body en- 
ergy level crossings which occur in the system at specific 
fillings and values of the magnetic fiux. The discontinu- 
ities of the persistent current seen in Fig. 1 correspond 
to the change of the total orbital momentum in the 2D 
array, which can be accompanied by spin-fiip transitions 
in the system. Finite disorder (W ^ 0) smoothens the 
discontinuities in the noninteracting (U, = 0) current, 
and removes some of the discontinuities in the interact- 
ing (U, V ^ 0) current. There are no discontinuities 
or spin fiips at half-filling in the interacting system (Fig. 
1(c)). (6) The sign of the persistent current at zero mag- 
netic fiux is generally dependent on the electron filling 
and the geometry of the array. We find the interacting 
(U ^ 0, V, W = 0) persistent current in 2 x 2, 3 x 2, 
4x2, and 3x3 arrays to be diamagnetic for n = 1/2_L 
and paramagnetic for n = 1/2. The random disorder of 
strength VF (= A, as used in our calculations) does not 
change the sign of the persistent current. (7) The current 
distribution in the clean 2D array is determined by the 
spatial symmetries of the system and band filling in the 
array. In the clean 3x3 lattice the current flows along the 
boundary for all electron flllings in the system. However, 
the persistent current distribution is disorder realization 
dependent. 

The most important feature of our results, the strong 
suppression of the persistent current at half-fllling by the 
on-site Coulomb repulsion [/, is a direct finite size mani- 
festation of the Mott-Hubbard metal-insulator transition, 
which is a true phase transition in the thermodynamic 
limit. The long-range Coulomb repulsion {V ^ 0) op- 
poses this transition in a flnite system (leading to an 
enhancement in the flnite U -suppressed current ampli- 
tude for y 7^ at half-fllling). Efl'ects of flnite V and W 
on the Mott-Hubbard transition in a 2D Hubbard model 
(in the thermodynamic limit) are not rigorously known. 



but our results are consistent with the expectation that, 
in general, flnite values of V and W should oppose the 
metal-insulator Mott transition at half-fllling. To further 
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FIG. 2. A plot of the RM^persistent current {Pfl"^ ver- 
sus the number of electrons N in the 3x3 array for the six 
different sets of U , V, W (as shown in the legends) in the 
3x3 system. All energies are given in the meV units. 

quantify the relative magnitudes of the effects of flnite 
U, V, W on the persistent current at various band flll- 
ings, we show in Fig. 2 our calculated root mean square 
current amplitude as a function of the number of elec- 
trons in the 3x3 quantum dot array for various values of 
Mott-Hubbard parameters. The rms current magnitude, 
^j2^i/2^ is obtained by averaging the typical current (Fig. 
1) over one flux period. The six different parameter sets 
with different values of U, V, W in Fig. 2 clearly show 
the dominant effect of the Mott-Hubbard transition (for 
flnite U) at half-fllling (N = 9), the generic destructive 
effect of flnite disorder at all flllings, and the generic ten- 
dency of the long-range Coulomb interaction homogenize 
the electron density in the system [12] and consequently 
to enhance the persistent current. 

In the well-studied single ID ring geometry, the effect 
of disorder (in the absence of any Coulomb interaction) 
on the persistent current is an exponentially strong sup- 
pression [13] of the current arising from the Anderson 
localization phenomenon associated with all electronic 
states being exponentially localized in one dimension in 
the presence of any flnite disorder. For our 2D quantum 
dot array the effect of flnite disorder (W ^ 0) is subtle 
(and much softer than in one dimension) because local- 
ization effects are logarithmically weak [14] in two dimen- 
sions for small W . (In fact, in the presence of an external 
magnetic held breaking the time reversal symmetry, the 
localization effect is even weaker [14] and may be hard 
to discern in a flnite sample.) This explains the rather 
benign effect of having a flnite W in Figs. 1 and 2 of our 
results, which should be contrasted with ID ring calcula- 
tions where the persistent current is strongly suppressed 
by flnite values of W . To better quantify our flnite disor- 
der results we show in Fig. 3 a log-log plot at half-fllling 
of the rms current {Py^^ without any interaction effects 
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FIG. 3. A log-log plot of {Pfl"^ versus W averaged 
over 100 disorder realizations for each value of W is shown 
([/ = y = 0) at half-filling in 3 x 3 (crosses), 4x4 (tri- 
angles), 5x5 (solid line), 6x6 (pluses). Asterisks show the 
interacting {U = 1 meV, = 0) results at half-filling for the 
3x3 system. Inset: {I^Y^^ at half-filling versus system sizes 
(i): triangles (2D array), squares (ID ring). (The average is 
taken over the range of the same total flux through the sys- 
tems, which varies from —4>ol2 to (l>o/2.) The dashed lines in 
the plot are the best linear fits to the corresponding data as 
described in the text. 

(U, V = 0), averaged over 100 disorder realizations for 
each value of W , as a function of the disorder strength 
W for various system sizes (3x3, 4x4, 5x5, 6x6). In 
plotting these results, we have factored out a scale factor, 
ric = (L^^^ — l^/L, so that the results for various sys- 
tem sizes fall on top of each other, showing approximate 
current scaling with system size and disorder. (The scale 
factor ric = (L^^^ — l)^ / L is the number of unit cells 
or plaquettes in each square array of size L.) All four 
non-interacting disordered results scale very well whereas 
the one disordered result with interaction deviates from 
the universal scaling at low disorder strength where the 
Mott metal-insulator transition effect dominates. The 
two dashed straight lines in Fig. 3 give the best fits to 
weak and strong disorder scaled currents, leading to the 
following empirical results for the effect of disorder on 
the persistent current: = (L/ric) (/(VF), with the 

scaling function g(W) being given by 

(W\ / 1^"^, 7 = (6-4±2.8) X 10-2, VF< 1.55 7rt 
'J(^>^ \W-l^, (]=1M±0A9, W>1.55Trt. 

The very small value of the scaling exponent j in the 
weak disorder limit is consistent with the expected loga- 
rithmic weak localization [14] in a 2D array, which even- 
tually crosses over to strong localization for large W. 
We also expect interaction effects to be particularly im- 
portant at weak disorder, but not so at high disorder, 
which is what is seen in Fig. 3 even at half-filling. Weak 
disorder, in fact, produces anti-localization effect (Fig. 
1(c) and 2) at half-filling m the presence of finite U , by 
slightly enhancing the persistent current from its finite U - 



suppressed value. Whether these last two findings remain 
true in the thermodynamic limit is, however, unknown. 

Note that the scale factor in Fig. 3 indicates that 
the persistent current scales with the perimeter size of 
the array or the square root of the number of dots, L^l"^ , 
in the system. We have explicitly verified this system size 
dependence by doing calculations for various system sizes 
at various fillings. As an example, we show as an inset 
of Fig. 3 a comparison of the system size dependence 
of the rms current in ID ring and 2D array geometries, 
finding an approximately linear dependence on L" with 
the best fit a 0.46 (1.1) in the 2D array (ID ring). 
Thus, the persistent current is proportional to the length 
of the boundary in both ID ring and 2D array geometry. 

In conclusion, we predict, based on an exact diagonal- 
izaton study of the generalized disordered Mott-Hubbard 
Hamiltonian, the existence of a ground state persistent 
current in finite 2D quantum dot arrays which should be 
experimentally observable. Direct observation of an equi- 
librium persistent current in coherent 2D quantum dot 
arrays will not only verify the formation of an artificial 
quantum dot molecule but will also shed light on the in- 
terplay among coherence, disorder, long- and short- range 
Coulomb interaction effects on Anderson-Mott-Hubbard 
quantum phase transitions. 
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